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Abstract 

A series of accurate ab initio calculations on CupOg finite clusters, properly embedded 
on the Madelung potential of the infinite lattice, have been performed in order to determine 
the local effective interactions in the Cu02 planes of La2_j;SrxCu04 compounds. The values 
of the first-neighbor interactions, magnetic coupling (JjvAr=125 meV) and hopping integral 
(tjVAf=-555 meV), have been confirmed. Important additional effects are evidenced, concern- 
ing essentially the second-neighbor hopping integral tjvAfJV=+110meV, the displacement of a 
singlet toward an adjacent colinear hole, h'gf)=-80 meV, a non-negligible hole-hole repulsion 
Vnn — Vnnn=0-Si eV and a strong anisotropic effect of the presence of an adjacent hole 
on the values of the first-neighbor interactions. The dependence of Jnn and tNN on the 
position of neighbor hole(s) has been rationalized from the two-band model and checked from 
a series of additional ab initio calculations. An extended t-J model Hamiltonian has been 
proposed on the basis of these results. It is argued that the here-proposed three-body effects 
may play a role in the charge/spin separation observed in these compounds, that is, in the 
formation and dynamic of stripes. 
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1 Introduction 



The insulating cuprates, such as La2Cu04, which are the undoped parent compounds of the high- 
Tc superconducting La2-xSr^Cu04, are known to present antiferromagnetic coupUngs between 
nearest neighbor (NN) copper centered sites in the Cu02 planes. Raman and neutron diffrac- 
tion experiments evaluate this coupling to be around J7VAr=130 meV(128ib6 meV[|^, 0,134ib5 
meV|^, 1^, ^, respectively). Nevertheless, the corresponding simple Heisenberg Hamiltonian does 
not reproduce entirely the features of the Raman spectrum ^, ^ |lO|, and additional 
effects such as second-neighbor magnetic coupling Jnnn and four-spin cyclic exchange have 
been invoked pi, P, ll|, |1. While an upper bond for Jnnn {{JnnnI ^9 meV) has been given 



from Raman experiments the amplitude of the four-spin operator in this kind of compounds 
is a matter of discussion. Previous works have shown that this cyclic operator corresponds to 
a fourth-order term in the Hubbard model, scaling as Xtj^j^/U'^, with A=40 [16, 17] or A=80 
[18, depending on the formal writing of the Hamiltonian. 

Regarding the hole-doped material, where the conduction takes place in the Cu02 planes, 
the holes can be seen as centered on copper atoms with large tails on the four neighboring oxygen 
atoms. They move from on site to an adjacent one through the effect of a hopping operator of 
amplitude tNN, for which there is no direct experimental evaluation, but values around -0.5 eV 



are considered as reasonable [20[. One of the most widely employed model Hamiltonians used 
for the interpretation of the properties of these materials, through a hole-pairing mechanism, is 
the so-called t — J model pl^, E2| which combines spin coupling and hole hopping: 



H = Y] Jnn • (Sa • Sb - 1/4) + tNN • {al^-b + alaa + s.c.)6{na + Ub, 1) (1) 



ab 



The adequacy of such a simple Hamiltonian to incorporate the physics of the problem is ques- 
tionable. Hopping between second-neighbor sites may be non negligible. When one derives the 
t — J Hamiltonian from the Hubbard Hamiltonian, three-site operators moving a singlet-coupled 
electron pair toward the hole appear at second order of perturbation theory, scaling as the Jnn 
operator, i.e, as tf^j^/U, where U is the on-site Coulombic repulsion. The transferability of 
Jnn from the undoped to the doped material is not guaranteed, the presence of a neighbor- 
ing hole may affect the coupling of two adjacent spins. The hole- hole repulsion Vij is likely to 
play a role, influencing the mean distance between the holes. Different extensions of the t — J 



model have been employed in numerical simulations, for instance, t — t' — J[23, 24, ^] or 
t — J — V^^, ^,but the values given to the parameters are rather arbitrary, varying widely 
from one author to the other, the main objective being to exhibit qualitative collective effects. 
Among them the experimental evidence of the occurrence of stripes have focused attention in 
the recent past [|30t[ . 



The goal of the present paper is to bring useful informations regarding the local effective 
interactions in undoped and hole-doped cuprates. To obtain them, the most accurate tools of 
ab initio quantum chemistry will be used. The method consists in considering few-site clusters, 
properly embedded in the field of the periodic environment, and to calculate the low part of the 
spectrum using the exact Hamiltonian, large basis sets and extensive Configuration Interactions. 
From this spectrum it is possible to fix the amplitudes of the effective interactions. The procedure 
has been successfully used to calculate the NN magnetic coupling in a series of perovskites[|3l]], 
among them La2Cu04 for which a value of JAr7v=130 meV is obtained. Similar calculations of 
the hopping integral in the hole-doped system (t=-0.57 eV) have been also reported []3^ , 
These calculations concerned symmetrical two-site clusters, for which the determination of Jnn 
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and tj\f]\f is straightforward from the two lowest eigenvalues. This is no longer the case when 
one considers larger clusters to extract additional parameters, concerning next-nearest neighbor 
(NNN) interactions, neighboring-hole dependence of Jnn and t^^, hole- hole repulsions and 
four-spin cyclic operators. Their achievement requires some additional mathematical tools such 
as the Bloch definition[^] of effective Hamiltonians and localization procedures (for instance. 
Boys method [^]). The methodology will be explained in Section II. Section III will present 
the results concerning the four-site square plaquette and three-site linear cluster with different 
number of electrons (holes). Calculations performed on undoped clusters provide the values of 
the NN, NNN and next NNN magnetic couplings (Scheme I). Also the amplitude of the four-spin 
cyclic exchange for this compound has been established. Hole-doped clusters give informations 
about the hopping integrals (NN, NNN and next NNN), the singlet-displacement operator and 
the dependence of the first-neighbor interactions (hopping integral and magnetic exchange) on 
the number and relative position of the adjacent holes. Section V presents a rationalization of 
the anisotropy of the effect of hole(s) in the vicinity on the values of Jnn and Inn and reports 
additional exploratory calculations to evaluate the dependence of the bicentric parameters on 
the hole positions. Finally, Section VI summarizes the results, proposing a refined t — J model, 
and discusses the possible effect of the additional operators on the charge/spin distribution on 
the lattice, with possible consequences for the stripping phenomena. 

SCHEME I 
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2 Method 

2.1 Mapping of a model Hamiltonian on an ab initio effective Hamiltonian 

For such materials the unpaired electrons are essentially located on Cu dx'^ — in-plane atomic 
orbitals, with non-negligible delocalization tails on the adjacent oxygen atoms. Such Cu-centered 
orbitals will be labelled {a, 5, c, ..}. In the doped material the hole has much larger delocalization 
tails on O 2p orbitals, but as shown by Zhang and Rice||2l|, it remains possible to work within 
a one-band model Hamiltonian, the precise nature of its valence orbitals being implicit. For a 
finite cluster involving p centers and n < p unpaired electrons, the model Hamiltonian works in a 
basis of n-electron localized determinants (pi, for which in the Heisenberg and t — J Hamiltonians, 
the double occupancy of the orbitals is prohibited. 

The ab initio calculations handle a large number of atomic orbitals, symmetry-adapted 
molecular orbitals (MOs) and expansions of the wave- functions on millions of determinants. 
Nevertheless, it is possible to construct from these calculations ab initio effective Hamiltonians 
which are in one-to-one correspondence with the model Hamiltonians. For an undoped cluster 
involving p Cu atoms {p sites), it is possible first to obtain from variational calculations a set 
of molecular orbitals containing doubly occupied MOs (core), unoccupied MOs (virtual MOs) 
and p MOs with essentially single occupation, which define the ab initio one-electron valence 
space, {(pi} in Figure 1. A unitary localizing transformation of the p symmetry-adapted orbitals 
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will provide equivalent localized orbitals {a',b',c'..} which can be seen as in strict one-to-one 
correspondence with the implicit valence orbitals of the model Hamiltonian. Fixing a double 
occupancy of the core orbitals, and putting n electrons in these orbitals, avoiding their double 
occupation, localized neutral determinants, {(j)'^ i^^}, are obtained which are in correspondence 
with the n-electron basis of the model Hamiltonian. These determinants define a model space (of 
projector P5, Ps = J2 \4'i ioc)i'f'i ioc\) initio calculations. Let be Ns the dimension of 

that space. The information obtained by the most refined ab initio calculations will be extracted 
according to the theory of the effective Hamiltonian proposed by Bloch |]3^ . When one knows 
the Ns eigenstates having the largest projections on the model space (which constitute the 
target space, stable subspace of Hexact) and their eigenvalues £'m, the effective Hamiltonian is 
such that its eigenvalues are the exact ones, and its eigenvectors are the projections of the exact 
eigenvectors onto the model space: 

H^f^\Ps^m)=Era\Ps'^ra). m=l,Ns (2) 

The spectral definition of W^^^ is: 

jjeff ^ ^ \Ps^>^)Em{Ps'^i\ (3) 

m 

where \Ps'^m) is the biorthogonal transformation of \Ps'^ni)- Actually the projections \Ps'^m) 
of the (orthogonal) states have no reason (expect for symmetry reasons) to be orthogonal, 
they define an overlap matrix s: 

Smn = {Ps'^m\Ps^n) (4) 

and the biorthogonal vectors are defined by 

\Ps^i)=s-^\Ps^^) (5) 

The values of the norms of the projections, i.e. the diagonal elements of s matrix, give an 
indication on the quality of the description of these states by the truncated space S. The model 
space and the exact eigenstates must be in strong correspondence, i.e. one must choose both 
spaces so that the vectors have the largest projections on the model space. 

Then one may express this Hamiltonian in the basis of the localized determinants {(p'^ 
written in terms of the orbitals {a', 6', c', .} and the matrix elements 

mH'^^W,) = Y.i'f>i\Ps^^n^)E^n^{Ps^i\c|>',) (6) 
m 

can be identified to the matrix elements {(j)i\H\(j)j) of the model Hamiltonian. In principle 
the effective Hamiltonians may be non-hermitian but the hermitization is straightforward 
The comparison between the ab initio effective Hamiltonian and the model Hamiltonian fixes the 
amplitudes of the integrals appearing in the latter and allows one to verify whether non-negligible 
additional interactions are not present. Figure 1 summarizes the whole process. Changing 
the size of the cluster, for instance going from a two-center/one-electron problem to a three- 
center /two-electron one, one may check the consistency of the procedure and the transferability 
of the effective interactions. 

2.2 Computational details 

The widely used embedded cluster technique has been employed to model the system. Finite 
clusters of the type CupOg (CU3O10 and CU4O12, Figure 2) have been selected, where the q 
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oxygen atoms are the first in-plane neighbors of the p Cu atoms. (Previous calculations have 
shown that the explicit involvement of the out of plane oxygen atoms does not change the 
values of the in-plane interactions ||3^, ^). The first-shell of atoms surrounding the cluster 
have been replaced by formal charges with pseudopotentials, in order to mimic the coulombic 
and exclusion effects. The rest of the lattice has been modeled by means of point charges, 
which values have been fixed according to Evjen's method|37|, and which correctly represent 
the Madelung potential of the crystal [ 38 1. This is a simplified approach in comparison with more 



elaborate methods proposed in the literature [ 39 1. 



The ten most internal electrons of Cu atoms have been represented by effective core pseudopo- 
tentials, the valence electrons being treated explicitly with triple-zeta basis sets. A double-zeta 
basis set has been used for oxygen atoms (preliminary calculations on bicentric clusters have 
shown that the inclusion of polarization functions on the bridging oxygen atoms has not an 
important effect neither on the magnetic coupling nor on the hopping integral) pO[| . 



For undoped clusters, the restricted-spin open-shell Hartree-Fock calculations variationaly 
define the singly occupied magnetic orbitals. These orbitals define a minimal valence complete 
active space (CAS). From this space it is possible to calculate the spectrum through a difference 
dedicated Configuration Interaction (DDCI) procedure which implies all the simple and 
double excitations on the top of this CAS, except the double excitations from the core to the 
virtual orbitals, which do not contribute to the energy difference at second-order of perturbation 



theory 



An alternative solution consists in defining an enlarged CAS including the on-bond 2p orbital 
of the bridging oxygen atoms. These ligand-centered orbitals are the most participating on the 



intersite spin-exchange and electron transfer processes [ 42 1 . Performing all the single excitations 
on the top of this extended CAS, which corresponds to the two-band Hubbard model, one 
introduces dynamical polarization effects, i.e., screening by the non-active electrons, at lower 
computational cost than the preceding computational scheme. 



3 Ab initio calculations on the plaquette and the linear clusters 

As was mentioned above, two different clusters have been used to extract the effective interac- 
tions. A four-site square cluster (plaquette) of formula CU4O12 has been employed in order to 
determine the first- (NN) and second-neighbor (NNN) interactions, and also the four-spin cyclic 
exchange. Third-neighbor interactions (nNNN) has been estimated by means of the calculations 
carried out in a linear three-site cluster (CU3O10). Comparing with the results obtained from 
the plaquette and previously studied binuclear clusters, it is possible to check the dependence 
of the NN interactions on the size of the fragment involved in the ab initio calculations. 



Three fillings of the valence shell have been considered in order to evaluate the dependence of 
these interactions on the hole concentration. Undoped (4centre/4electron and 3centre/3electron 
problems), one hole-doped (4c/3e and 3c/2e) and two hole-doped (4c/2e and 3c/ le) situations 
have been analyzed. From the systems with two holes in the valence shells, it is possible to 
extract the amplitude of the hole-hole repulsions, an important magnitude for the study of 
the hole pairing mechanism. It is worth to notice that the here-referred hole-dopings are not in 
correspondence with the total doping of the lattice, induced by the replacement of La"^^ by Sr"*"^. 
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A change in the occupation of the valence shell of these small clusters just induces a local hole 
doping, which provides informations about the local modifications of the effective parameters. 



3.1 The localization process 

In the plaquette, the four symmetry-adapted valence orbitals belong to the Aig, and Big 
representations in the D^^^ group. The localizing unitary transformation is straightforward since: 

aig = {a + h + c + d)/2 (7) 

e«(i) = {a + h-c-d)/2 (8) 

eu{2) = {a-b-c + d)/2 (9) 

big = (a-b + c-d)/2 (10) 

Figure 3 pictures one of these four localized valence orbitals for undoped and for the doped 
plaquettes, showing the strong localization of the magnetic orbitals and the d — p hybridization 
occurring in the hole-doped systems] 32, 33, 43]. 



In the linear cluster the localizing transformation of the three magnetic orbitals is no longer 
imposed by the symmetry. Two of these orbitals {ipg^ip'g) belong to the Ag irreducible repre- 
sentation and the other one {(pu) to the Biu symmetry. The rotation [/, which transforms the 
{ifg^ip'g^ipu} into the localized set {a', ?)',(/} has been performed according to the Boys crite- 
rion |35], which maximizes the distance between the centroids of the orbitals. An alternative 



localization criterion, the minimization of the direct exchange integral Kac-, leads to the same 
rotation. 



3.2 The magnetic interactions 

Table 1 reports the results obtained from the plaquette and the linear clusters for the magnetic 
coupling, involving different fillings of the valence shell. For undoped cluster, the first-neighbor 
magnetic coupling Jnn is ~ 125 meV. The calculated value is independent on the size of the 
considered cluster and it is in agreement with both the previously determined Jnn in binuclear 
cluster and the estimations from Raman and neutron diffraction experiments (128ib6 meV[j|, 
^, 134ib5 meV[P, ^ ^, respectively). The coupling between second-neighbors Jnnn is also 
antiferromagnetic, with a value of Jnnn=Q-5 meV in accordance with the experimental upper 
limit \ Jnnn\ <9 meV [^. A negligible antiferromagnetic coupling has been found between Cu 
atoms placed at a 2R distance (third-neighbors, next NNN): JnNNN=^ meV. 

Finally, the calculations on the plaquette provide a non-negligible value for the four-spin 
cyclic exchange of K= 14 meV. The four-body operator can be written as: 

Hk = K [{Si■Sj){Sk■Sl) + {Si■Sl){Sj■Sk)-{S^■Sk){SJ■Sl)] (11) 

<ijkl> 

This operator produces the cyclic permutation of the four spins on the plaquette (Scheme II) 
plus ordinary two-spin exchanges of all the pairs of spins of the plaquette including those of 
the diagonals. (A detailed discussion has been reported elsewhere p^j.) Its value is somewhat 
smaller than some estimations used in numerical simulations []l4| , |13| , |45|l , but larger than the 
critical value, (ivT/ J^v at )c=0. 05 ±0.04, estimated by Sakai and Hasegawa [15| for the appearance 
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of a magnetization plateau at half the saturation value in the S = ^ antiferromagnetic spin 
ladders. 

SCHEME II 

^1 \l — *l 

t' \ 

When a hole is introduced in the cluster, the NN magnetic coupling is influenced by its 
presence, but in different directions depending on the relative position of this hole. Thus a 
hole in a colinear position to the two NN spins increases the coupling between these two spins 
(Scheme Ilia). However the NN magnetic coupling diminishes if the hole is placed in a position 
perpendicular to the bond (Scheme Illb). The same trend is observed in the plaquette when a 
second-hole is introduced, Jnn being 94 meV, to be compared with Jnn= 104 meV in presence 
of 1 hole and J7VAr=125 meV for undoped systems. 
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3.3 Hopping integrals 



As in the case of magnetic interactions, the first-neighbor hopping integral {t^N) is also inde- 
pendent on the size of the cluster (Table 2), a value of -558 meV has been found in Ihole-doped 
linear cluster, -552 meV in Ihole-doped plaquette and -555 meV in previously reported binuclear 
cluster This value is in accordance with a generally accepted value of -500 meV for these 



compounds 1 20 1 



An evaluation of second- and third-neighbor hopping integrals has also been possible. The 
NNN hopping integral t]sf]\fN=+112 meV is unexpectedly large. The sign is in agreement with 
the negative overlap of the active orbitals placed at a \/2R distance, but its magnitude is large 
due to through bond processes, which involve the oxygen atoms. On the basis of the perturbation 
theory, there are two contributions to the t^NN hopping integral. One corresponds to a a third- 
order contribution, scaling as where t,, is the hopping integral between the O 2p 

and the Cu 3d orbitals, tpp is the hopping integral between the O 2p orbitals and AEct is the 
O 2p — > Cu 3d charge-transfer excitation energy (Scheme IV). Since tpp, tpd and AEct are 
negative quantities, the third-order contribution results in a positive magnitude. There exists 

an alternative pathway corresponding to a fourth-order contribution, scaling as .^i , with 



CT 



opposite sign. Since the third-order contribution is expected to be larger than the fourth-order 
one, the resulting sign of t^NN is positive. 

A moderate amplitude of the hopping integral between third-neighbors (distance 2R) has 
been extracted from the calculations on the linear cluster. The extension of this coupling is 
controlled by through-bond interactions, but in contrast with the plaquette, only fourth-order 
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process can be written here, smaller in magnitude and with an opposite sign with respect to 

tNNN- 

SCHEME IV 
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When an additional hole is introduced in the system, the NN hopping integral is unchanged 
when the hole is placed in a perpendicular position to the bond (t,/vjv=-558 meV in the 2hole- 
doped plaquette, Scheme Va), but its absolute value is augmented when the hole is placed in a 
colinear position to the bond (tArAr=-600 meV in the 2hole-doped linear cluster. Scheme Vb). 
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The presence of the extra hole does not influence significantly neither the second-neighbor 
itNNN= +130 meV to be compared with +112 meV in absence of this additional hole) nor the 
third-neighbor {tnNNN=-'i^ meV versus -47 meV) hopping integrals. 



3.4 Singlet displacement operator 

An additional information coming from these calculations concerns the singlet-displacement 
operator(Table 3). It is a three-site/two-electron operator, which moves the pair of electron, 

coupled in a singlet, toward a hole placed in a neighbor position. Thus, a singlet on sites a and 
b, c containing a hole, is displaced to the positions b and c, the hole being in a (Scheme VI): 

SCHEME VI 




It can be written as: 

hfr) [\ab - ab){bc - bc\ + \bc - bc){ab - ab\] 6{na + m + Uc, 2) (12) 

where 6{na + rif, + nc, 2) controls the fact that the three centers bear only two electrons, i.e. 
that the singlet can only move toward a hole. As in the precedent parameters, the amplitude 
hg'j^ depends on the relative position of the hole. If the singlet moves to a neighbor bond in the 
plaquette (i.e., the hole is placed in a perpendicular position to the singlet bond), the value is 
/i|'g= -41 meV while it goes to /i|'g=-80 meV if the hole is on the same axis than the singlet. 

The presence of a second hole in the plaquette does not affect the amplitude of the singlet- 
displacement {h'g£,= -37 meV). 

There exists also a small similar operator involving the diagonal of the square 

hf^ [1 6c - bc){bd- bd\ + \bd-bd){bc - bc\] (5(nf, + + n^, 2) (13) 

moving the electrons as shown in Scheme VII. The value of this operator in presence of a second 
hole is very close to the preceding value: hg'^= 14 meV. 

SCHEME VII 




3.5 Hole-hole repulsions 



The absolute value of the hole-hole repulsion is not accessible, but the effective Hamiltonian 
gives the difference between two situations (Table 4). Prom the plaquette, we can extract the 
relative stability of two holes placed in NNN positions with respect to two holes in NN: 

Vnn - Vnnn = Vab - Vac = 0-98eV (14) 



This is significantly larger than the values usually accepted for simulations [27, 28, |29[]. One 
should stress on the fact that this value takes into account the dynamical repolarization effects 
of all the atoms explicitly treated in the calculation, i.e., the screening by the twelve in-plane 
oxygen atoms linked to the four Cu atoms of the plaquette. It misses the polarization of the 
rest of the environment. Taking into account the polarization of a large surrounding shell, with 
a value of the polarizability of the ion, ao = 1.30 (obtained from a series of finite-field 
ab initio calculations on an embedded CuO^ cluster), in agreement with experimental estimates 



|46|, one diminishes the hole-hole repulsion difference to Vnn — VAr7VAf=0.80 eV which remains 
a rather large value. From the linear cluster, it is possible to estimate the energy gain obtained 
when placing two holes in nNNN positions (distance 2R) instead of adjacent positions (distance 
R): Vnn — VuNNN = Vab — Vac= 1-77 eV, with a final value of Vnn — VnNNN='^-^^ eV once the 
environment polarization effects have been taken into account. 

These values of differences between hole-hole repulsions may seem very large and most of 
the calculations introducing this repulsion in a t-J-V model use to take smaller values (J < 



9 



V < 4J)[p^, 28, 29| when trying to exhibit hole-pairing phenomenon. The above calculated 



values are smaller than the corresponding electrostatic quantities calculated in a point charge 
approximation: Vatat — VArAr7v=l-12 eV and Vnn — VnNNN=^-^0 eV. The delocalization of the 
holes on the oxygen atoms should result in larger repulsions, especially for Vnn since the two 
holes share an oxygen atom, and our values, which exhibit a significant screening, do not seem 
unrealistic. 



4 Interpretation of the hole dependence of the one bond Jjy^ 
and tNN parameters 

The two preceding sections indicate that the presence of holes in the immediate vicinity of a 
bond may affect the values of the spin coupling and of the hopping integral on this bond, and 
that this effect is anisotropic, i.e., does not only depend on the minimal distance of the hole to 
the atoms of the bond. Hence it seems important to rationalize these effects if possible. 



4.1 Rationalization from the two-band model 



The hopping integral. 

The rationalization of this anisotropy is possible in terms of the two-band model. The t hopping 
integral results from a second order effect (Scheme VIII, where the open circle represents a 
bridge oxygen atom and closed circles correspond to copper atoms): 

SCHEME VIII 
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t - (15) 

where AEqt is the 2p O to 3d Cu charge transfer excitation energy. In presence of an additional 
hole, the energies of the model space determinants, (pi and 4'2, and of the intermediate charge 
transfer state, (pcT, will be modified. The model space determinant energies are no longer 
degenerate and one shall take their mean energy as the zero-order energy. Let consider now 
the effect on the bond directed along x of an adjacent hole placed either on the y direction, 
as occurs in the plaquette, or on the x axis, as in the linear cluster (schemes Va and Vb). In 
situation (Va) the zero-order mean electrostatic energy with the hole is 6Eo = |^=^ = In 
the corresponding charge transfer state the mean electrostatic energy with the hole is SEqt = 
~ excitation energy is increased by a small quantity: 

l^E'cT = AEcT + SEcT - SEq = AEcT + (16) 

which should diminish slightly the absolute value of the hopping integral. In contrast, for 
situation (Vb), 6Eo = ^ and SEct = ^ hence A£^^-p = AEqt — which should increase 
the absolute value of t, in agreement with our calculation. 
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This analysis is quite rudimentary, it neglects the possible effects of the polarization of the 
orbitals on the tpd integrals and rests on a crude evaluation of the electrostatic effects on the 
denominators (holes considered as centered on Cu sites, neglecting their oxygen character), but 
it seems to agree with the computed trends. 

The magnetic coupling. 

The same kind of analysis may be tempted for the Jnn coupling constant, which is now a 
fourth-order quantity (Scheme IX): 

•^^^ ~ UiAEcTy ^^^^ 

which involves intermediate situations of opposite polarities, but with the same energy: AEqti = 
AEcT2 = AEcT- 111 the undoped system, Ui=U2=U, which corresponds to the Coulombic re- 
pulsion of two electrons placed in the same 3d orbital. 

SCHEME IX 
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Adding an external hole will not modify the electrostatic zero-order energy, but will affect the 
energies of all intermediate states, so that the coupling becomes: 

^^^'^ ■ (for the top pathway) + ^^j"' • ^ (for the bottom one). 

Notice that while AEcti = AEqt + SEcti and AEct2 = AEqt + 5Ect2 have no reason to be 
related, Ui = U + 6U and U2 = U — SU whatever the outer charge distribution. The effective 
coupling is therefore: 

_,4 ( 1 1 I 1 (18) 

P'^ V {AEcT + SEcTif U + SU {AEcT + SEcT2f U - SU J ^ ^ 

To the first order in ^ and developments, one gets: 

•^^^ - UiAEcT? \ " AE^ " AE^ ~ W ^ W ) ^ ^^^^ ^^^^ 

The effect of the additional holc(s) will go through their electrostatic interaction in the charge 
transfer states. For the perpendicular situation (Scheme X) 6Ecti = and 5Ect2 = 
Hence, dEcTi + SEct2 = and Jnn is diminished by the presence of the adjacent hole in 
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the y direction, as observed in the corresponding ab initio calculations (AJjvAr/«^ArAr=-17%). 
SCHEME X 
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Oppositely, in the linear situation (Scheme XI) the quantities are: SEcti = ^ and SEct2 = 
so SEcTi + SEcT2 = ^ is a negative quantity, hence the coupling constant should be in- 
creased. This is actually observed since AJjVAr/Jjv7v=+25%. The variation of Jnn in the model 
is of the right sign and its amplitude is larger than for the perpendicular orientation. 

SCHEME XI 
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T/ie singlet- displacement operator. 

The same kind of rationalization applies to the two-clectron/three-centre operator h'g'^. In the 
one-band model this effect scales as t'^/U' (Scheme XII): 
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SCHEME XII . . 
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where U' = U — 2Vab + Vac, Vij being the electrostatic interaction between the charges (with 
respect to the basic undoped distribution). The denominator will be smaller for the linear con- 
figuration than for the perpendicular one since Vac is larger (l/\/2ii instead of 1/2R). This may 
explain the difference between the amplitudes of Ngj^ found in the plaquette and in the linear 
cluster (Table III). The more realistic two-band model introduces an alternative mechanism, 
since there exist two paths to go from ab to be, cf Scheme XIII: 
SCHEME XIII 



t tl i t I t I t| t i 

• y • o • « » . • o • e • 

+ + - + 

^ AEi AE2 



^ ^ AE3 

t Mi H \\ ti \\ 

• o • e • — *- • o • e • 

+ - + + - + 

AE'i AE'2 
The second path corresponds to the process appearing in the one-band model, which is not the 
case for the first one. Explicitly, the excitation energies are: 



A£;i = Ed-Ep-Up = AEcT-Ud + -^ (20) 



2 
R 

AE2 = 2Ed-2Ep + Up-2Up-^ + Vi,n = 2AEcT-Ud + Vi,ii (21) 

AE3 = Ed-Ep-Up + Ud-^ + Va,ii = AEcT-^ + Va,ii (22) 
AE[ = AE3 (23) 

^E'2 = Ud-^+Vac (24) 

It is likely that AE2 > AE2, i.e. that the additional path has a larger contribution. The overall 
effect is: 

iprf / 1 1.1 1 



AE[ \AEi AE2 ^ AE[ ' AE'^J ^^^^ 

ie, of the same order of magnitude as Jnn- Regarding the orientation effect it is clear that the 
denominators are larger for the perpendicular orientation than for the linear one, due to the 
hole-hole repulsions Viji,Viji and V13, and it is expected that hg^^^ will be larger for the linear 
orientation, as found in the numerical ab initio calculations (-80 meV and -40 meV, respectively). 
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4.2 Further exploration of the influence of the hole on J^n and t^^ 

In order to evaluate this influence, larger clusters would have to be considered, which is beyond 
our computational posibilities (the preceding CI expansions frequently reach 10^ determinants), 
and a simpler procedure has to be used. In the present set of calculations, a dinuclear CU2O7 
clusters have been chosen as before, adding one or two additional point charges on Cu sites of 
the environment, which become Cu^^ centers, in order to grossly mimic the effect of the holes. 
For first-neighbor holes, this is a rather crude approximation, and its relevance has to be tested 
by comparison with the previously reported four- and three-center calculations. Actually, for 
a perpendicular position of the hole JArAr=115 meV (compared to 104 meV in the plaquette). 
For a colinear hole Ji\f]\[=175 meV, compared to 156 meV in the linear cluster. Hence the 
anisotropy of the first-neighbor hole is qualitatively reproduced. As an additional confirmation 
we have calculated the effect of two holes in the plaquette and we get the same value JArAr=93 
meV as for the 2hole-containing plaquette. Table 5 gives the results for Jnn concerning a series 
of possibilities with one remote hole, two holes in different positions, together with the variation 
of the electrostatic energy in the intermediate charge transfer states with respect to the non- 
doped case, expressed in eV. Figure 4a shows the correlation between the calculated Jnn value 
with the variation of the electrostatic energy SEqti + SEct2 in the intermediate charge transfer 
states. The correlation is quite satisfactory and suggests the following law in eV: 

Jnn = 0.131 - 0m3{5EcTi + SEct2) (26) 

An analogous fit of the hopping integral dependence on the variation of the electrostatic energy 
in the intermediate charge transfer state 5Ect — ^Eq has been tempted from the values of tjsiN 
(Table 5) calculated on binuclear complexes in presence of external hole(s). The correlation is 
less satisfactory (Figure 4b). The distortions of the active orbitals in presence of these close holes 
should be non-negligible and should affect the amplitudes of the hopping integral in a rather 
complex manner. However, it seems simpler to research a linear law rather than to produce and 
handle a dictionary of operator amplitudes considering exhaustively all possible occurences. We 
therefore propose the following fit in eV for the first-neighbor hopping integral: 

= -0.521 + 0.187(<5^cT - SEo) (27) 



5 Discussion and Conclusions 

The numerical results obtained here above suggest that the usual t-J or t-J-V model Hamil- 
tonins neglect important physical effects, which should be incorporated into an extended model 
Hamiltonian. The main deviations from the t-J Hamiltonian are the following: 
(i) the inclusion of the second-neighbor hopping, t^NN — +110 meV. The third-neighbor in- 
teraction tnNNN — -40 meV could tentatively be omitted. Recent analysis of angle-resolved 
photoemission spectroscopy data has shown that both tjyNN and tnNNN are necessary for un- 
derstanding the dispersion and line shape of the spectral function in the t-J model pH]. Tohyama 



et al. [25| have estimated the ratio tj^NN /tNN and tnNNN /tNN to be -0.12 and 0.08, respectively, 
by fitting the tight-binding Fermi surface to the experimental one in the overdoped sample [^] 
on the assumption that in the overdoped region the Fermi surface of the tight-binding band is 
the same as that of the t-t'-t"-J model. The values here-proposed are in good agreement with 
these ratio: t^NN /tNN =-0.19 and tnNNN /tNN=0-08 but not with the ratio J/t^N proposed 

by these authors {J /tNNTohyama=^-'^ J/tNN=0.22). 
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(ii) the singlet displacement operator hsD has to be taken into account, at least for the colinear 

displacement, since hsD — -80 meV. 

(iii) the hole-hole repulsion appears to be far from negligible. Due to possible screening effects 
and electrostatic mean cancellations by the Sr ions, it is certainly reasonable to neglect the hole- 
hole repulsions beyond the third-neighbors. Even if the calculated values Vnn — Vnnn=^-8 
eV and Vnn — ^nJVAriv=l-5 eV were somewhat exagerated, these interactions certainly play an 
important role. 

(iv) the magnetic coupling and hopping integral between adjacent atoms depend in a stere- 
ospecific manner on the existence and position of hole(s) in the immediate vicinity. A simple 
correlation with the electrostatic energies of the intermediate ligand to metal charge transfer 
states has been proposed resulting in simple formulas, which should be used in a realistic model 
Hamiltonian. Prom the exploratory calculations appearing in Table 5, and for sake of simplicity, 
it seems sufficient to consider the first neighbors of the bond in the calculation of the electro- 
static energy changes SEcti + SEct2 and SEct — ^Eq appearing in formulas (27) and (28), 
respectively. 

One may therefore propose the following extended t-J model Hamiltonian: 

NN 

H = *(<^^) ■ (^1^1 + "-W + s-c-)S{np + Uq, 1) - J{SE) ■ {Sp ■ Sq) + 

<pq> 
NNN 

+ 51 ^JVATAT • (flpflr- + + S.C.)S{np + Ur, 1) - JnNN ' {Sp ■ Sr) + 

<pr> 

plaquette 

+ £ ■ [(5p ■ Sq) {Sr ■ Ss) + {Sp ■Ss){Sq-Sr)- {Sp ■ Sr){Sq ■ S,)] + 

<pqrs> 
connected 

+ ^ hFg^ ■ (a^ataqOf + a^a^qaqap — a^pol^araq — aJatopOg + s.c.)5{np + Uq + rir, 2) -|- 

<pqr> 
<nNNN 

+ J2 Vpq-5{np,0)-6{nq,0) 

pq 

In this Hamiltonian, t{SE) and J{SE) reflect the dependence of the first-neighbor interactions 
on the number and position of the adjacent holes, which is controlled by the energy changes in 
the charge transfer intermediates. The second-neighbor interactions are not influenced by the 
presence of adjacent holes and the displacement of a singlet takes placed toward an adjacent 
hole, so the positions occupied by the singlet and the hole have to be connected. 

The here-proposed modifications with respect to the usual t-J or t-J-V Hamiltonians are 
important. We would like to point out that some local physical effects evidenced in the present 
work may have an impact on the spatial ordering of charges an spins and on their dynamics. In 
a crude static look at this problem one may notice that stripping separation of charges and spins 
(i) favors the mobility of holes in the charged column, since t and J values become especially 
larger, (ii) leads to a stronger magnetic coupling in the bonds perpendicular to the stripes ( with 
a trend to form singlets on these bonds). The two last effects would result from the three body 
corrections on t and J created by the adjacent holes. Finally, (iii) the mobility of the stripes can 
be enhanced by the second-neighbor hopping integral and by the singlet displacement operator. 

We think that our calculations are reliable cnoTigh and that their qualitative conclusions and 
quantitative estimates deserve to be considered in simulations of the collective properties of the 
lattice, whatever the method used for such a study. 
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6 Figure captions 



Figure 1. Summary of the strategy used to extract effective interactions in La2-a;Sr2;Cu04 sys- 
tems. 

Figure 2. Representation of the clusters employed in the ab initio calculations, (a) the plaquette 

CU4O12, (b) the linear cluster CusOio- The first-neighbors of the cluster atoms have been also 
included. They have been modeled by means of the combination of a pseudopotential and a 
point charge, to mimic the both the exclusion and the Coulombic effects. 

Figure 3. Localized active orbitals for the plaquette: (a) for undoped systems, centered in 
3dx2-y2, with important tails on the four in-plane neighbor oxygen atoms, (b) for hole-doped 
systems a strong 3d — 2p rchybridization takes places, the 2p character increases sustantially 
with respect to the undoped situation. 

Figure 4. Linear correlation between the calculated first-neighbor interactions in binuclear clus- 
ters and the change of the electrostatic energy in the intermediate charge transfer states (in eV) : 
(a) Jnn versus SEcti + SEct2', (b) Inn versus 6Ect — SEq. 



19 



Table 1. Magnetic interactions in undoped and hole-doped clusters (in meV). 



Undoped clusters 



Exp. 



Doped clusters 



J. 



NN 



t i 
I 1 



• »-^^ — I 



122 



I t 



124 



128±6*, 134±5'= 



• t-i — • 



156 



I + 



104 



A K 



+ + 



94 



Jnnn 



I ■■••■■0 



6.5 



< 9° 



CK + 



10 



J, 



nNNN 



-11 



K 



'■■•4.. 



14 



a. A value of Jjvjv=130 meV has been previously reported [32,33], where polarization functions have 
been included in the basis set of the bridge oxygen atom. b. Ref. [1,2]. c. Ref. [3,4,5]. d. Ref. [5]. 
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Table 2. Hopping integrals. Dependence on the presence of an additional hole (in meV). 



tNN 



tNNN 



1'nNNN 



IHole- doped clusters 



-555" 



-558 





-552 



+112 



+ 



-47 



2Hole- doped clusters 



-600 




-558 



+130 



-36 



a. Ref. [32,33]. 
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Table 3. Singlet-displacement operator (in meV). Effect of the presence of a second hole in 
the neighborhood of the singlet. 



1 Hole- doped clusters 



''SD 




-41 



''SB 



i V. 

• •- 



-80 



+ 
-• 



2 Hole- doped clusters 




+ + 



-37 
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Table 4. Hole- hole repulsions (in eV). 



Vnn — Vnnn 


Vnn - VnNNN 




i + + 1 t + 


1 1 1 - 1 1 1 

1.77(1.47) 




+ 


+ 






1 • 


1 + + i 

0.98(0.8) 
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Table 5. Magnetic coupling and hopping integral in binuclear clusters with hole(s) in adjacent 
position(s) and intermediate state energies (all in eV). 





Undoped 


IHole-doped 


J 


SEcTl + ^EcT2 


t 


SEcT ~ SEq 












0.116 


0.305 


-0.450 


0.152 
























0.135 


-0.17 


-0.558 


-0.084 




















0.175 


-0.61 


-0.590 


-0.305 




















0.093 


0.61 


-0.458 


0.312 




























0.098 


0.61 


-0.482 


0.305 


























0.116 


0.57 


-0.456 


0.305 


























0.130 


-0.14 


-0.502 


0.076 






















0.144 


-0.685 


-0.593 


-0.335 
























0.193 


-1.22 


-0.691 


-0.633 
























0.200 


-0.84 


-0.563 


-0.404 
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Figure 1. Calzado and Malrieu 
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Figure 2. Calzado and Malrieu 



26 



0.25 




0.05 ' — ' ' — ' ' ' — ' ' — ' ' ' 

-1.5 -1 -0.5 0.5 1 

5E,,,+5E,,,(eV) 




Figure 4. Calzado and Malrieu 



28 



